N PS  ARCHIVE 
1966 
KEITH,  W. 


■11 


J9QBBS 

MlMufltlimM 


IHHN 


A  STATISTICAL  MODEL  FOR  DETERMINATION  OF 
THE  RADIATIVE  TEMPERATURE  AT  "BLACK- 
BODY'  SURFACES 


WILLIAM  H.  KEITH 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unlimit< 


. 


A  STATISTICAL  MODEL  FOR  DETERMINATION  OF 
THE  RADIATIVE  TEMPERATURE  AT  "BLACK-BODY" 
SURFACES,  BASED  UPON  SPECIFIC  WINDOW- 
RADIANCE  VALUES  MEASURED  OVER  A  RANGE  OF 
ZENITH  ANGLES  BY  NIMBUS  II  RADIOMETER 

by 

William  H.  Keith 
Lieutenant  Commander,  United  States  Navy 
B.A.,  University  of  Wisconsin,  1953 

Submitted  in  partial  fulfillment 
for  the  degree  of 

MASTER  OF  SCIENCE  IN  METEOROLOGY 

from  the 

UNITED  STATES  NAVAL  POSTGRADUATE  SCHOOL 
October  1966 


MPs  A?ci\\\JC 


ABSTRACT 

Computed  Channel  2  or  window  radiances  for  the  Nimbus  II 
optical-sensor  system  were  kindly  made  available  to  the  author  by 
the  Meteorological  Satellite  System  of  ESSA  for  each  of  the  106 
model  atmospheres  considered  by  Wark  et  al.  (1963) 

Multiple  regression  equations  relating  the  emitted  black-body 
intensity  from  "black"  interfaces  to  the  independent  variables  (1) 
"observed"  specific  filtered  radiance  and  (2)  total  optical  path, 
were  set  up  for  each  of  seven  zenith  angles  in  the  range  ©=0°  to 
60°,  and  for  each  of  62  randomly  selected  cases  from  the  Wark  sounding- 
catalog.  Both  variables  gave  high  statistical  significance,  with  the 
former  accounting  for  the  primary  part  of  the  variance  of  the  dependent 
variable,  but  with  variable  (2)  always  contributing  in  such  a  way  as 
to  account  for  some  of  the  atmospheric  absorption  along  the  sensing 
di  recti  on. . 
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1.  Introduction 

Wark  et  al.  \j~^   ,  and  other  authors  have  shown  that  both  the 
surface  air  temperature  under  clear  sky  conditions  and  the  cloud-top 
temperature  when  the  sky  is  covered  by  relatively  dense  clouds,  as 
stratocumulus  or  altostratus,  can  be  estimated  using  outgoing  terres- 
trial radiation  intensity  observations  obtained  by  Tiros  II,  III  and 
IV  channel  2  radiometers.  Channel  2  radiometers  are  filtered  to 
receive  only  infrared  radiation  in  the  "window"  region  of  the 
spectrum,  and  in  the  Nimbus  II  this  window  region  was  10-11/4  . 
Usually  the  estimated  temperature  was  several  degrees  too  low  under 
clear  sky  conditions. 

By  utilizing  Nimbus  II  channel  2  radiometer  computations  of 
sensed  outgoing  terrestrial  radiation  intensity  (as  computed  by  the 
Meteorological  Satellite  Laboratory)  for  62  atmospheric  models  -- 
all  of  the  47  cloud-undercast  cases?  and  15  randomly  selected  clear- 
sky  cases,  multiple  linear  regression  equations  are  developed  to 
solve  for  the  angular  dependence  of  total  outgoing  specific  intensity 
1(9  )  values.  The  multiple  linear  regression  analysis  uses  the 
BMD03R  program  as  adapted  to  the  CDC  1604  computer. 

In  order  to  describe  radiative  temperatures  at  the  earth's 
surface  under  clear  sky  conditions,  and  at  the  cloud-top  for  under- 
cast  conditions  an  estimate  of  the  total  specific  intensity  as  viewed 
by  Nimbus  II  must  be  made  as  a  function  of  the  emitted  black  body 
intensity  and  of  an  absorptivity  dependent  upon  the  effective  optical 
path  {L^SEC  0     ).  Since,  in  general,  window  radiation  is  more 
readily  available,  it  was  decided  to  use  specific  filtered  (window) 
radiation,  rather  than  total  specific  intensity.  The  latter  may  then 


be  inferred  as  a  function  of  interface  temperature  and  optical  path, 
as  is  done  in  Section  7. 

Graphs/tables  are  then  constructed  which  allow  simple  compu- 
tation of  the  surface  air  temperature  under  clear  skies,  and  cloud 
top  temperatures  using  total  water  vapor  mass,  zenith  angle,  and 
satellite  radiometer  measurements  in  Channel  2. 

2.  Data  Sources 

The  data  utilized  consisted  of  the  outgoing  terrestrial  radiation 
intensity  computed  by  the  Meteorological  Satellite  Laboratory  from 
window  region  (Channel  2)  data  obtained  from  the  Nimbus  II  satellite. 
106  atmospheric  models  were  developed  by  the  Meteorological  Satellite 
Laboratory  in  order  to  calculate  the  outgoing  intensities  (_7J  ,  of 
which  62  selected  cases  are  used  in  this  paper.  The  62  cases  are 
composed  of  47  cases  undercast  at  various  altitudes  and  15  clear 
cases,  which  are  identified  in  Appendix  I.  The  models  developed 
by  the  Meteorological  Satellite  Laboratory  were  taken  from  radio- 
sonde data  over  all  latitudes  and  seasons.  All  the  soundings  went 
to  at  least  25  mb,  and  the  temperatures  were  extrapolated  parallel 
to  a  standard  atmosphere  to  the  pressure  of  O.lmb.  Stratospheric 
humidity  was  extrapolated  to  conform  in  general  to  standard-atmosphere 
values.  The  sky  condition,  the  temperature  at  selected  levels  and 
the  temperature  of  the  undercast  level,  as  well  as  the  total  water 
vapor  and  ozone  masses  were  also  listed  as  part  of  the  sounding 
catalog  in  \JJ  . 

In  order  to  deduce  total  emergent  flux  and  the  angular  dependence 
of  these  flux  estimates,  the  specific  filtered  intensity  in  the  10-11/4 


band  for  each  of  the  62  atmospheric  models  was  used  at  each  of  the 
four  given  zenith  angles  (0=0°,  20°,  45°  and  60°).  The  inten- 
sities at  three  additional  zenith  angles  {0   =  10°,  32.5°  and  52.5°) 
were  obtained  by  interpolation. 

The  radiation  intensities  of  the  62  models  were  given  in  ERG/ 
(CM^  SEC  STRDN)  and  are  converted  to  watts  per  meter  per  strdri 
in  this  thesis. 

3.  Nature  of  the  Data,  and  of  the  Model 

The  Nimbus  II  satellite  uses  both  a  High  Resolution  Infrared 
Radiometer  (HRIR)  scanning  the  4.0  micron  window,  and  a  Medium 
Infrared  Resolution  Radiometer  (MRIR)  which  scans  the  infrared 
spectrum  in  regions  where  the  far  infrared  region  of  the  solar 
spectrum  gives  only  insignificant  contributions.  Only  Channel 
2  radiation  of  the  MRIR  segment  is  used,  and  the  sensed  radiation 
comes  from  an  area  of  approximately  40  miles  across  in  the  hori- 
zontal (based  upon  a  half-cone  viewing  angle  of  2.9°).  The  solid 
angle  thus  presented  at  the  sensor  is  sufficiently  small  so  that  each 
measurement  must  be  considered  a  specific  intensity  rather  than  flux. 

The  emergent  flux  is  a  more  important  quantity  than  that  of 
outgoing  specific  intensity  from  the  meteorological  point  of  view, 
because  the  latter  applies  to  a  specific  direction  Q  only,  whereas 
the  emergent  flux  gives  total  radiative  heat  loss  from  a  column. 
In  atmospheric  radiation  problems  it  is  convenient  to  treat  the 
ground  and  cloud  surfaces  as  black  bodies  \J\j   .  As  shown  in 
Q3  pp.  81-82J  the  intensity  and  flux  of  radiation  from  a  black  body 


are  respectively 

in  terms  of  the  black-body  surface  temperature  Tg.  In  (1),  the 
constant  0~  of  the  Stefan-Boltzman  law  has  the  value 

a- =0.5 6  6  -y  X  /07watts/m* 

so  that  all  values  of  FR  are  in  watts  /rrr  and  those  of  Ig  are  in 
watts /m^/steradi an . 

The  emergent  filtered  intensity,  I  ^  (©),  scans  primarily 
the  10-11/M  window  region.  The  nature  of  the  filter-transmissivity 
function  in  this  band  is  an  important  consideration.  The  techniques 
by  which  the  filter- functions  are  employed  in  conjunction  with  the 
sensor  absorptivity  is  discussed  in  considerable  detail  by  Wark 
et  al  [7,  p.  3,  and  Supplement,  p.  2_J.  These  same  considerations 
have  been  applied  to  the  measured  values  of  I  ^  {Q)   for  Nimbus  II, 
and  the  computations  provided  by  the  Meteorological  Satellite  Lab- 
oratory are  presumably  available  to  other  interested  workers.  The 
interesting  facts  are  that  (1)  the  thermister  sensitivity  is  not 
unity  in  any  of  the  wave  band  intervals,  (2)  that  window  filter 
function  is  not  unity  over  the  window,  nor  is  it  zero  outside  this 
spectral  region.  However,  suffice  it  to  say  Wark's  filtered  radiance 
(_7,  Eq.  3 J  multiplied  by  the  effective  thermistor  emissivity 
provides  the  "observed  window  radiation"  of  this  study.  Such  data 
is  of  the  form 


10 — // 

5 


iv  — . — a 


(10 


with  C  known  as  a  quasi-constant  with  respect  to  I  q  {&)   over  a 
wide  angular  range.  Theoretically  one  could  then  enter  a  table  of 
*  Q)x&)  ~   va^ues»  divide  by  C,  and  recover  the  appropriate  temper- 
ature for  (jj~h   )  from  the  Planckian  monochromatic  intensity  function. 

W  / 

However,  even  in  the  so  called  window  region  there  are  small  atmospheric 
depletions  due  to  both  water  vapor  and  ozone,  which  should  be  restored 
if  accurate  "black-surface"  temperatures  are  to  be  estimated.  Thus, 
in  (3)  below  a  depletive  effect  is  added  back  to  recover  part  of 
the  initial  beam  intensity. 

Hence  the  linear  regression  was  set  up  with  a  heuristically 
physical  basis,  as  follows: 

77*  u        © 

Here  M*  has  been  multiplied  by  a  single  temperature-ratio  correction 
factor  of  the  Lorentz-type,  since  the  individual  computations  of  U* 
did  not  include  such  factors  in  the  layer-by-layer  summations.  The 
inclusion  of  this  temperature-ratio  term^  along  with  the  constant  A2» 
(to  be  determined)  was  aimed  at  inducing  a  more  homogeneous  strati- 
fication of  the  extinction  treatment,  since  the  cloud  top  cases 
investigated  here  varied  considerably  with  respect  to  pressure  level. 
Thus,  (3)  becomes 

q^y=  A>+A,I» (&)+AtA  (I)  I/e)  i±S£C 6    (4» 

Moreover  in  view  of  the  one-to-one  (expected)  relationship  between 
I  ^  (  Q)   and  ("TT^  remarked  upon  above,  a  linear  statistical 
relationship  of  form 
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Tje)*J*,(*-r')+A,t 


(5) 


is  proposed,  with  both  k-j  and  l<2  quite  small.  In  fact,  apart  from 
sensor-noise,  and  second-order  depletive  effects,  l<2  should  be 
negligibly  small  and  is  treated  thus  here.   (4)  then  assumes  the 
form 

where  A0,  A-j ,  A  ,  and  k  would  each  have  to  be  determined  by  standard 
best-fit  methods.  The  grouping  Tnl  $   (  G  )\  is  emphasized  in  order 
to  show  the  conversion  of  intensity  to  flux  units  throughout  (6). 
We  therefore  arrive  at  an  alternative  form  of  (6) 

Y, '«  b,  v-  SeV-^Xs  (7) 

where  the  independent  variables  are 

x3=<r-r'(x/\sece 

and  the  dependent  variable  is 

Y,  =^7"¥  (8b) 

The  aim  of  the  working  equation,  (7),  is  to  determine  how  well 
the  effective  black-body  temperature  can  be  described  statistically 
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(8a) 


by  use  of  observed  values  of  filtered  radiance  Iaj(  ©)  and  of  a 
more  complex  variable  X  ,  which  appears  to  depend  upon  the  dependent 
variable,  but  actually  its  variance,  Q"x^  £w  CON$T[  O"   jC  ^  J  $ 
However  Xg  was  formed  as  shown  with  values  of  U.  also  known  in  advance 
for  each  of  the  62  soundings  used,  in  accordance  with 

(9) 


P>*  0,1*8 


The  choice  of  the  standard  temperature  T  =  300°K  was  used  in 
connection  with  the  full  Lorentz-type  broadening  factor  as  set  back 
into  (3).  This  value  was  adopted  based  primarily  upon  experimental 
transmissivities  of  Palmer  ~5J  who  found  that  he  could  reproduce 
the  theoretical  absorption  coefficients  of  Yamamotc  as  listed  in[_7j 
computed  for  T  =  260°K.  This  difference  is  well  brought  out  by  Wark 
et  al .  I  7,  eq.  5j  ,  which  suggests  that  in  order  to  retain  U*  in  the 
general  format  of  (9),  the  exponent  of  (p/p0)  should  be  slightly 
less  than  unity.  But  then  if  strict  conformity  with  Lorentz-type 
broadening  is  required,  the  standard  temperature  T0  is  raised  to 
300K. 

4.  Data  Processing 

As  data  samples  for  the  testing  of  (7)  at  seven  different  angles 

in  the  range  Q=   0°  to  60°,  the  cloud-top  and/or  surface  temperatures 

4 
were  used  as  input  for  Y-j  =  0~T  for  each  of  the  62  soundings  chosen 

from  M.S.L.  Report  Number  10  [_7J  .  For  each  such  sounding,  values 

of  I  q  (  9  )  at  angles  6  =  0°,  20°,  45°,  60°,  were  also  given  for 
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Nimbus  II.  The  I q   (  &  )  sample  was  arbitrarily  enlarged  by  linear 
interpolation  to  corresponding  values  at  10°,  32.5°, 52.5, giving 
seven  angles  for  each  sounding  with  seven  corresponding  values  1^  (  &  ). 
The  variable 


x3=(^vXiooKst-ce 


was  formed  by  transgeneration  using  a  program  available  for  this 
purpose  at  U.S.  Naval  Postgraduate  School  Computer  Facility.  There 
were  7  collections  of  punched  data  cards  for  each  sample  set  (Y2.X2.X3), 
each  collection  corresponding  to  a  different  zenith  angle  applied  to 
the  set  of  62  soundings. 

The  multiple  regression  program  BMD03R,  as  adapted  to  the 
CDC  1604,  was  utilized  to  solve  for  the  value  of  the  best  fit 
coefficients  for  each  of  the  seven  zenith  angles.  The  total  input 
to  the  program  consisted  of  7x62  or  434  data  cards.  Table  1  shows 
a  sample  input  for  one  zenith  angle. 

5.  Statistical  Results  and  Inferences 

The  multiple  linear  regression  equations  obtained  from  the  program 
are  shown  in  Table  2,  in  the  form  of  (7)  and  (8).  It  can  be  readily 
seen  by  examining  the  formulas  of  Table  2,  that  the  regression  coeffi- 
cient of  the  first  independent  variable  considerably  over-weights 
that  variable  when  compared  to  the  regression  coefficient  of  the 
second  variable,  the  depletion  factor. 

In  order  to  examine  the  question  of  a  significant  specification 

of  CTT4  by  the  independent  variables  X0  and  X3=  CTT^i)**  u  sec  8  , 

c  *       lo 
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Table  1 

Sample-set  Employed  in  the 
Multiple  Regression  Analysis  * 


e 

X2  =  i  ©  ( e  ) 

X3(watts/M2) 

Y1(watts/M2) 

10 

17.116 

40.084 

262.854 

10 

22.008 

108.177 

323.953 

10 

17.672 

76.523 

275.150 

10 

26.423 

495.751 

373.705 

10 

13.087 

11.423 

220.381 

10 

15.787 

40.668 

247.104 

10 

27.962       1071.279 

406.182 

10 

16.987 

50.299 

262.854 

10 

24.897 

399.156 

363.327 

10 

21.474 

129.534 

319.265 

10 

25.129 

181.637 

363.327 

10 

8.635 

1.694 

158.511 

*The  above  12  samples  represent  12  of  62  samples  for  the  zenith 
angle  9  =  10°. 
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Table  2 


Multiple  Linear  Regression  Equations  Obtained  from 
the  BMD03R  Program  for  Seven  Zenith  Angles 


e 


i 


Dp    Xp 


B3  X3 


10u 

20° 

32.5° 

45° 

52.5° 

60° 


srrJ 
ctt4 

TT4 

O-T4 
O-T4 


58.270 

+ 

11.9226X2 

+ 

58.869 

+ 

11.8840X2 

+ 

52.412 

+ 

12.3788X2 

+ 

62.265 

+ 

11.9233X2 

+ 

56.224 

+ 

12.2479X2 

+ 

57.978 

+ 

12.1866X2 

+ 

63.657 

+ 

12.1443X0 

+ 

01001X. 


01084X. 


00724X. 


00998X. 


00983X. 


01125X- 


01121X. 
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the  F-statistics  in  Table  3  will  be  used  as  obtained  from  the  output 
of  the  multiple  linear  regression  program. 

The  BIMD03R  program  tests  for  significance  of  variables 
introduced  into  a  multivariate  linear  regression  equation.  The  most 
readily  adaptable  test,  in  the  analysis  of  variance,  is  the  so-called 
F-test  where  the 


F-statistic  =  Mean  squares  explained  by  predictor      (iq) 
Mean  squares  unexplained 


In  this  experiment  X~  is  the  primary  specifying  variable,  and  for 

its  test  F  =  F,  ^  is  defined  as 
1 ,60 

p    =  Sum  of  squares  explained 
1 ,60  fSum  of  squares  unexplained! 

L        wrr  J 

The  subscript  numbering  (1,60)  denoted  the  number  of  degrees  of 
freedom  associated  with  numerator  and  denominator,  respectively. 

The  addition  of  the  depletion  factor,  symbolized  by  X3,  requires 
the  definition  of  a  new  F-statistic  specifically  restricted  to  the 
testing  of  X3,  after  the  contribution  of  X^  has  been  excluded.  This 
statistic  will  be  denoted  F  and  be  defined  as 

p'    _  Additional  Sum  of  squares  explained         (n) 
1,59   fSum  of  squares  unexplained"] 
L        (N-3)  J 

The  critical  value  at  the  99%  confidence  level  of  belief  is  7.08  for 
both  F  and  F  ,  and  these  critical  values  are  surpassed  at  the  99% 
confidence  level  for  each  variable  added,  with  the  exception  of 
X3  at  B=   20°  and  32.5°.  However  the  ratio 


F'.*o(**l  =  0(V)  to    0  do3) 
7.06 
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TABLE  3 

Some  Statistics  Relating  to  the  Multiple  Linear 
Regression  Equations  for  Seven  Zenith  Angles 


Partial         Mult  Computed 

Reg.Coeff  Corr  Coeff  Corr  Coeff  F  &  F'   Values 


9     Vari-     Mean 
able 


19.690 
272.446 
296.471 


11.92261  .99669 

.01001  .55242 


9987 


23927.718 
25.912 


10( 


20( 


Yl 

x3 
Yl 


32.5°  X? 
*3 
Yl 

45°    X0 


Tl 
52.5°  X2 


60° 


19.681 
275.661 
296.471 

19.620 
289.959 
296.471 

19.436 
323.062 
296.471 

19.248 
385.321 
296.471 

19.131 
458.220 
296.471 

18.720 
544.308 
296.471 


11.88413    .99463 
.01084    .49824 


12.37883 
.00724 


11.92329 
.00998 


12.24789 
.00983 


12.18661 
.01125 


.96464 
.14609 


.96001 
.24503 


99519 
61318 


,98884 
56546 


12.14434    .94758 
.01121    .36280 


9980 


9859 


9876 


9982 


9959 


9798 


14765.000 
18.482 


2088.534 
1.287 


2972.724 
3.769 


22457.562 
35.550 


7231.164 
27.733 


1427.700 
8.943 
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and  is  so  far  beyond  critical  that  the  contribution  of  X~  is  highly 
significant,  and  certainly  cannot  be  accounted  for  by  a  chance 
contribution  in  any  of  the  seven  cases.  The  contribution  of  X3 

astestb*FWx^\   i 

7708    ' 
in  all  cases  except  &  =  20°  and  32.5°. 

It  should  be  noted  that  the  coefficient  of  the  depletion 
factor  X  is  positive  in  all  cases;  that  is,  it  helps  to  account 
for  a  small  but  significant  part  of  absorption  from  the  originating 
flux  CTT  ,  which  part  has  been  absorbed  in  the  atmosphere.  It  has 
therefore  been  tentatively  concluded  that  all  7  angular  regression 
formulas  afford  significant  information  with  respect  to  the  atmo- 
spheric depletion,  as  well  as  to  the  primary  variable  I  (p  (  &  ). 
Here,  in  effect  some  weight  is  being  given  to  the  sign  test  of  the 
regression  coefficient  of  X_,  even  for£  =  20°  and  32.5°,  where  the 
F1  test  falls  short.  The  F-test  failure  in  the  case  &  =   32.5  is  due 
to  the  inclusion  of  Q  =  20°  data  in  the  formation  of  the  &    =  32.5 
sample. 

An  explanation  for  the  lack  of  significance  of  the  depletion 
factor  f or  0  =  20°  could  lie  in  the  interplay  of  water  vapor  and 
ozone  radiations,  weak  as  both  of  these  are  likely  to  be. 

The  shortcoming  of  the  statistical  model  must  be  recognized: 
only  one  absorbing  element  has  been  introduced,  namely  water  vapor, 
which  absorbs  weakly  in  the  "window"  region)  7,  fig.  3J  of  the 
atmosphere.  In  addition,  the  small  ozone  content  introduces  a 
moderate  absorbing  band  where  the  water  vapor  concentration  is 
sparse  (15  km  and  above).  The  emergent  intensity  formula  for  a 
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single-absorber  atmosphere  is  given  by 

T,(e)'T8A^h^?(^ec^(AT)     (12) 

and  in  a  wave  number  interval  where  two  absorbers  give  overlapping, 
the  transmissivity  "fcrj  is  generally  taken  as  the  product  of  the 
individual  transmission  functions  for  the  gases  taken  separately. 

First  of  all  it  appears  that  in  the  direct  computation  of  top- 
transmitted  radiance,  there  may  be  an  angle  where  overlap-effects 
give  larger  discrepancies  than  at  extremely  close  or  extremely 
distant  angular  ranges.  Secondly,  the  statistical  attempt  of 
synthesizing  two  absorbers  by  putting  in  just  water  vapor  could 
not  account  for  all  of  the  variance  attributable  to  depletion. 
Minor  variations  in  the  proportionality  constant  C  of  (2)  undoubtedly 
do  occur  for  any  given  setting  of  filtering  wave  band  limits. 

On  the  whole,  however,  a  physically-based  statistical  model 
which  has  a  multiple  correlation  coefficient  »  0.98  is  quite 
adequate  for  specification  of  the  dependent  variable,  and  in  the 
present  case  of  accounting  for  approximately  one- third  of  the 
unexplained  variance  left  over  when  Ia  {&)   alone  is  used  in 
specification. 

6.  Estimation  of  Surface  and  Cloud  Top  Temperatures 

Utilizing  the  multiple  regression  equations  listed  in  Table 
2,  for  the  seven  zenith  angles,  Fig.  1  and  Appendix  II  have  been 
developed  to  afford  close  approximations  of  the  surface  and/or 
cloud-top  temperatures.  The  basic  input  data  for  use  in  these 
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graphical  aids  are  the  Nimbus  II  Channel  2  radiation  measurements, 
together  with  the  water-vapor  path  in  direction  &   .     These  graphs 
are  based  upon  (15)  which  is  derived  below. 
Since  by  (7) 


VT^^  B.-i-Bz**^3*** 


we  may  define,  a  unique  black-body  flux  using  the  mean  values, 


B,  and  B- 


(TT^fli^X, 


5 


(13) 


Then  making  the  approximation^"^  =  J*  Tp  in  X3,  we  have 

q.rf.0-Ar+-i/'-a*^«*\*.  — 11  (14) 


5      3 


f(<ri»<&)V«»J  " 


By  use  of  the  binomial  expansion  carried  out  to  the  linear  term  only 
the  water-vapor  corrected  radiative  temperature  is  therefore  given 
accurately  by 


r=v 


From  (15),  a  first  estimate  of  the  surface  and/or  cloud  top 
temperature  is  obtained  based  upon  the  averaged  form  of  the  regression 
equation  (7).  This  is  justifiable  insofar  as  (13)  is  concerned  since 
the  coefficients B,  and  B2  vary  from  their  means  by  only  5%,  an  insig- 
nificant deviation  from  the  mean. 
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The  bracketed  factor  in  (15)  may  be  written  as  a  depletive 
correction  to  the  interface  temperature  ATR  (  Q  ) 

&(&)-%  z2*L)  %£ce 


vn 


t 


Values  forAT„(  ©  )  are  listed  in  Tables  4,  5 10,  for  6  =0°, 

D 

10 ......  60°,  respectively.  The  input  data  for  each  table 

is   ■  &(Qj\       and  the  optical  depth  14*. 

Tests  of  the  tables/graphs  against  known  interface  tempera- 
tures show  wery   good  results.  The  corrections  A  TB  are  less  than 
2K,  and  usually  less  than  0.5K,  but  always  additive  in  the  sense 
of  adding  back  depleted  flux.  In  fact,  the  two  steps  just  described, 
namely  forming 

%+ATB 

results  in  specification  of  the  interface  temperature  T  accurate  to 
+  0.2K.     An  alternative  method  for  obtaining  ATB  is  discussed  in 
Section  7. 


7.     Angular-Dependence  Corrections  and  Conclusions 

The  major  results  derived  have  been  the  determination  of  a 
composite,  pooled, first  estimate  of      TB  from  (13).     Since    Tg     has 
been  derived  by  averaging  all   of  the  nearly  identical   forms  of  the 
terms  B-j+B2X2  of  regression  equation  (7),     TB  may  be  considered 
independent  of  the  zenith  angle. 

The  depletive  correction  AT${@/  to  the  radiative  temperature, 
arising  from  the  presence  of  the  water-vapor  term  was  of  the  form 

A-£(e)  =  4(e)^<4 
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in  the  regression  equation  (7,  15).  Here  k  (Q  )  is  a  constant 
for  each  selection  of  ©  ,  but  differs  from  one  zenith  angle  to 
the  next. 

To  emphasize  the  nature  of  the  angular  dependency  which  now 
resides  solely  in  the  depletion  term  &lg(&)   ,  consider  the 
cases  9   =   0°  and  60°:  ^ 

9-0°       ATB(Uir)o)  =  Z,3b(>TsQ>X/d 
©=60°     AT8(utfSEC(>6)~5'.Z52l%Ut*<6:i 

This  specialization  is  done  not  only  because  these  results  span 
completely  the  angular  range  employed,  but  also  because  both  cases 
corresponded  to  significant  F'-tests  for  variables  X  in  the  analysis 
of  variance  (Section  4). 

For  intermediate  values  of  Sj^'sl"/   may  be  estimated  by 
the  Taylor  expansion  formula  applied  in  the  form 

atb(&)  =  ^T5^)+[|2^]a(s£c©)  '      (i6) 

and  when  finite-differences  are  employed, 


xrB  Ce) = atb  (o) +r^B(eo)-^Ta(o)lper  e  _ ,) 

rom  which  we  finally  obtain 


(17)  fits  the  results  of  Appendix  II,  Tables  4  -  10,  within  +  0.2K 
and  may  be  considered  a  good  estimator  of  the  best-fit  depletive- 
correction  for  all  .angles  9  in  the  range  0°^  ©  ^    60°. 
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Further,  the  combination  of  (17)  together  with   TD  of  (13)  led  to 

D 

the  recovery  of  the  known  interface  temperature  to  an  accuracy  of 
+  0.2K. 

As  a  measure  of  the  effects  of  energy  depletion,  consider 
the  intensity-depletion  /±H{Q) 


or 


With  regard  to  flux  absorbed  in  the  atmosphere,  relative  to  the 
cloud-top  estimate,  FB(  T  +^1T),  one  may  integrated  I  (_    ©  ) 
hemi spherically  to  obtain      ^ 

o 

Upon  integration,  the  last  equation  gives  the  flux-loss  between 
interface  and  the  channel   2  sensor 

(19) 


^f  =  Fb(^b)^  \jOSl  X  /0"3 


Here  the  last  expression  in  the  bracket  shows  the  result  of  the 
integration.  Limb-darkening  has  been  neglected,  since  according  to 
Wark  et  al  |  1  \    only  about  25%  of  the  possible  flux  lies  outside 
of  a  60°  angle  cone. 

Even  though  the  cloud  temperature  sensed  by  Channel  2  may  be 
too  low  by  0.5°  to  2.0°,  depending  upon  the  zenith  angle,  the 
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"sensed"  black-body  intensity  and  flux  values  at  cloud  tops  are  too 
low  by  approximately  1  to  2%   depending  upon  the  actual  value  of 
/yjZ\B)/To      of  (17)  and  upon  the  value  U*  in  (19),  respectively, 

Refinements  of  this  statistical  model  may  be  made  to  include 
a  posteriori  values  of  ozone  or  of  dust-scattering,  but  such 
refinements  are  not  likely  to  be  as  statistically  significant  on 
an  operational  basis  where  the  optical  masses  of  ozone  and  of  dust 
are  not  known  until  long  after  the  data  readout  is  made. 
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APPENDIX  I 
15  Clear  Sky  Atmospheric  Models  from  [__7J 
Used  in  Multiple  Linear  Regression  Program 


Number 

1 

2 

3 

4 

7 

8 
10 
12 
13 
20 
23 
27 
31 
50 
51 


Location 
Havana,  Cuba 
Kenitra,  Morocco 
Oakland,  California 
Nantucket,  Massachusetts 
Ft.  Smith,  Arkansas 
Thule,  Greenland 
Chaquaromos  Bay,  Trinidad 
Tucson,  Arizona 
Washington,  D.C. 
Clark  Field,  P.I. 
Stephenville,  Texas 
Glascow,  Montana 
Barter  Island,  Alaska 
Chateauroux,  France 
Bitteurg,  Germany 


Time 

1200  GMT,  Sept.  29,  1958 

1200  GMT,  Sept.  29,  1958 

1200  GMT,  Sept.  29,  1958 

1200  GMT,  Sept.  29,  1958 

1200  GMT,  Sept.  29,  1958 

1200  GMT,  Sept.  29,  1958 

1200  GMT,  April  1,  1958 

1200  GMT,  April  1,  1958 

1200  GMT,  April  1,  1958 

1200  GMT,  Jan.  1,  1958 

1200  GMT,  April  1,  1958 

1200  GMT,  Jan.  1,  1958 

1200  GMT,  Jan.  1,  1958 

1200  GMT,  Jan.  25,  1958 

1200  GMT,  Jan.  25,  1958 
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APPENDIX  II 
CORRECTIONS  TO  FIRST  ESTIMATES  OF  TB  BASED 
ON  WATER  VAPOR  PATH  FOR  ZENITH  ANGLES0=1O°, 
32. 5°,  45°,  52.5°,  and  60°  (Tables  4,  5,  6, 

7,  8,  9  and  10) 
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TABLE  4 


Corrections  to  First  Estimates  of  Tg  Based 
on  Water  Vapor  Path  for  Zenith  Angle 


e  =  o' 


^ 

0.1 

0.2 

0.4 

0.8 

1.6 

3.2 

0 

.0424 

.0848 

.1696 

.3392 

.6784 

1.3568 

2 

.0462 

.0924 

.1848 

.3696 

.7392 

1.4784 

4 

.0493 

.0986 

.1972 

.3944 

.7888 

1.5776 

6 

.0517 

.1034 

.2068 

.4126 

.8272 

1.6544 

8 

.0540 

.1080 

.2160 

.4320 

.8640 

1.7280 

10 

.0561 

.1122 

.2244 

.4488 

.8976 

1.7952 

12 

.0579 

.1158 

.2316 

.4632 

.9264 

1.8528 

14 

.0595 

.1190 

.2380 

.4760 

.9520 

1.9040 

16 

.0611 

.1222 

.2444 

.4888 

.9776 

1.9552 

18 

.0620 

.1240 

.2480 

.4960 

.9920 

1.9840 

20 

.0638 

.1276 

.2552 

.5104 

1.0208 

2.0416 

22 

.0650 

.1300 

.2600 

.5200 

1.0400 

2.0800 

24 

.0662 

.1324 

.2648 

.5296 

1.0592 

2.1184 

26 

.0673 

.1346 

.2692 

.5384 

1.0768 

2.1536 

28 

.0684 

.1368 

.2736 

.5472 

1.0944 

2.1888 

30 

.0695 

.1390 

.2780 

.5560 

1.1120 

2.2240 

32 

.0704 

.1408 

.2816 

.5632 

1.1264 

2.2528 

34 

.0712 

.1424 

.2848 

.5696 

1.1392 

2.3784 

36 

.0722 

.1444 

.2888 

.5776 

1.1552 

2.3104 

38 

.0731 

.1462 

.2924 

.5848 

1.1696 

2.3392 

40 

.0739 

.1478 

.2956 

.5912 

1.1824 

2.3648 

30 


TABLE  5 

Corrections  to  First  Estimates  of  TR  Based 
on  Water  Vapor  Path  for  Zenith  Angle 


e=  10 


2^ 

0.1 

0.2 

0.4 

0.8 

1.6 

3.2 

0 

.0467 

.0933 

.1866 

.3732 

.7464 

1.4928 

2 

.0508 

.1016 

.2032 

.4064 

.8128 

1.6256 

4 

.0541 

.1082 

.2164 

.4328 

.8656 

1.7312 

6 

.0569 

.1138 

.2276 

.4552 

.9104 

1.8208 

8 

.0594 

.1188 

.2376 

.4752 

.9504 

1.9008 

10 

.0615 

.1230 

.2460 

.4920 

.9840 

1.9680 

12 

.0635 

.1270 

.2540 

.5080 

1.0160 

2.0320 

14 

.0653 

.1306 

.2612 

.5224 

1.0448 

2.0896 

16 

.0670 

.1340 

.2680 

.5360 

1.0720 

2.1440 

18 

.0685 

.1370 

.2740 

.5480 

1.0960 

2.1920 

20 

.0800 

.1400 

.2800 

.5600 

1.1200 

2.2400 

22 

.0713 

.1426 

.2852 

.5704 

1.1408 

2.2816 

24 

.0726 

.1452 

.2904 

.5808 

1.1616 

2.3232 

26 

.0738 

.1476 

.2952 

.5904 

1.1808 

2.3616 

28 

.0750 

.1500 

.3000 

.6000 

1.2000 

2.4000 

30 

.0760 

.1520 

.3040 

.6080 

1.2160 

2.4320 

32 

.0871 

.1542 

.3084 

.6168 

1.2336 

2.4672 

34 

.0782 

.1564 

.3128 

.6256 

1.2512 

2.5024 

36 

.0792 

.1584 

.3168 

.6336 

1.2672 

2.5344 

38 

.0802 

.1604 

.3208 

.6416 

1.2832 

2.5664 

40 

.0811 

.1622 

.3244 

.6488 

1.2976 

2.5952 

31 


TABLE  6 


/\ 


Corrections  to  First  Estimates  of  T     Based 
on  Water  Vapor  Path  for  Zenith  Angle 

o=zo° 


0.2 


0.4 


0.8 


1.6 


3.2 


0 

.0314 

.0628 

.1256 

.2512 

.5024 

1.0048 

2 

.0346 

.0692 

.1384 

.2768 

.5536 

1.1072 

4 

.0371 

.0742 

.1484 

.2968 

.5936 

1.1872 

6 

.0391 

.0782 

.1564 

.3128 

.6256 

1.2512 

8 

.0410 

.0820 

.1640 

.3280 

.6560 

1.3120 

10 

.0425 

.0850 

.1700 

.3400 

.6800 

1.3600 

12 

.0440 

.0880 

.1760 

.3520 

.7040 

1.4080 

14 

.0452 

.0904 

.1808 

.3616 

.7232 

1.4464 

16 

.0464 

.0928 

.1856 

.3712 

.7424 

1.4848 

18 

.0476 

.0952 

.1904 

.3808 

.7616 

1.5232 

20 

.0486 

.0972 

.1944 

.3888 

.7776 

1.5552 

22 

.0496 

.0992 

.1984 

.3968 

.7936 

1.5872 

24 

.0505 

.1010 

.2020 

.4040 

.8080 

1.4160 

26 

.0513 

.1026 

.2052 

.4104 

.8208 

1.6416 

28 

.0521 

.1042 

.2084 

.4168 

.8336 

1.6672 

30 

.0530 

.1060 

.2120 

.4240 

.8480 

1.6960 

32 

.0536 

.1072 

.2144 

.4288 

.8576 

1.7152 

34 

.0545 

.1090 

.2180 

.4360 

.8720 

1.7440 

36 

.0551 

.1102 

.2204 

.4408 

.8826 

1.7652 

38 

.0557 

.1114 

.2228 

.4456 

.8912 

1.7824 

40 

.0563 

.1126 

.2252 

.4504 

.9008 

1.8016 

32 


TABLE  7 


/\ 


Corrections  to  First  Estimates  of  TR  Based 
on  Water  Vapor  Path  for  Zenith  Angle 


0  =  3ZS 


im^ 

0.1 

0.2 

0.4 

0.8 

1.6 

3.2 

0 

.0508 

.1016 

.2032 

.4064 

.8128 

1.6256 

.2 

.0550 

.1100 

.2200 

.4400 

.8800 

1.7600 

4 

.0586 

.1172 

.2344 

.4688 

.9376 

1.8752 

6 

.0615 

.1230 

.2460 

.4920 

.9840 

1.9680 

8 

.0641 

.1282 

.2564 

.5128 

1.0256 

2.0512 

10 

.0663 

.1326 

.2652 

.5304 

1.0608 

2.1216 

12 

.0685 

.1370 

.2740 

.5480 

1.0960 

2.1920 

14 

.0704 

.1408 

.2816 

.5632 

1.1264 

2.2528 

16 

.0722 

.1444 

.2888 

.5776 

1.1552 

2.3104 

18 

.0737 

.1474 

.2948 

.5896 

1.1792 

2.3584 

20 

.0753 

.1506 

.3012 

.6024 

1.2048 

2.4096 

22 

.0768 

.1536 

.3072 

.6144 

1.2288 

2.4576 

24 

.0781 

.1562 

.3124 

.6248 

1.2496 

2.4992 

26 

.0793 

.1586 

.3172 

.6344 

1.2688 

2.5376 

28 

.0807 

.1614 

.3228 

.6456 

1.2912 

2.5824 

30 

.0818 

.1636 

.3272 

.6544 

1.3088 

2.6276 

32 

.0829 

.1658 

.3316 

.6632 

1.3264 

2.6528 

34 

.0841 

.1682 

.3366 

.6732 

1.3464 

2.6928 

36 

.0852 

.1704 

.3408 

.6816 

1.3632 

2.7264 

38 

.0862 

.1724 

.3448 

.6896 
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